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Abstract 

We consider a particle in the over-damped regime at zero temperature under the influence of a sawtooth 
potential and of a noisy force, which is correlated in time. A current occurs, even if the mean of the noisy 
force vanishes. We calculate the stationary probability distribution and the stationary current. We discuss, 
how these items depend on the characteristic parameters of the underlying stochastic process. A formal 
expansion of the current around the white-noise limit not always gives the correct asymptotic behaviour. 
We improve the expansion for some simple but representative cases. 
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1 Introduction 



In the last two years there has been a considerable interest in the problem of noise induced transport. 
The motivation to study such models has been initiated by Magnasco [j|]. He showed that there are 
two necessary ingredients for noise induced transport: A stochastic force that is correlated in time and 
an environment without inversion symmetry. In the one-dimensional case one can have e.g. a periodic 
potential without inversion symmetry. Such a potential is often called a ratchet-like potential. 

Magnasco was mainly interested in the adiabatic limit, i.e. in the case where the fluctuations of the 
stochastic force are slow. A few month later, the problem was investigated more carefully by Doering et 
al. [g]. They studied the motion of a one-dimensional particle in a saw-tooth potential and in the case 
where the force is described by various stochastic processes. For the symmetric dichotomous Markov 
process and at vanishing temperature they obtained an exact solution for small correlation times. For 
the Ornstein-Uhlenbeck process and for a class of processes called kangaroo processes, they calculated 
the current to first order in the correlation time. They also presented some results from Monte-Carlo 
simulations of such a system, again at zero temperature. Later, Mielke [Q] developed a method that 
allows to calculate the current for a large class of processes including the ones discussed by Doering et 
al. [jUl, again in the case of a saw-tooth potential or more generally for a piecewise linear potential. He 
recovered the results from []3|] and found several other cases where a current reversal occurs. One of 
these cases was a process that consists of an even sum of dichotomous processes. For such a process his 
results differ from the perturbative result to first order in the correlation time [ph. The sign of the current 
was different, even for very small values of the correlation time. One of the motivations of the present 
work is this discrepancy. We will show that the current and the stationary distribution of the coordinate 
of the moving particle are non-analytic functions of the correlation time. To do this we generalize the 
method developed in [|J] so that it is applicable to zero temperature as well. This allows us to obtain exact 
numerical results for the current and the stationary distribution, as well as a correct asymptotic expansion 
for small correlation times. This is also of general interest, since there are many approximative methods 
to solve Langevin or Fokker-Planck equations in this or other cases that agree with the usual perturbative 
expansion for small correlation times (for a review we refer to Furthermore we show that some 
of the properties of the solution for small correlation times are relevant for larger correlation times as 
well. Even if the solution is analytic for small correlation times, it may become non-analytic for larger 
correlation times. This always happens if the support of the stationary distribution of the stochastic force 
is finite. 

The paper is organized as follows: In the next section we define the model and the stochastic processes 
we will be able to deal with. In the subsequent section we show how the method in [Q] can be generalized 
so that one can treat the zero temperature case as well. Since a large part of the method is similar to the 
finite temperature case, we refer to [Q] for details. In section 4 we present some of the numerical results 
including a detailed discussion of the current reversal for various noise processes. In section 5 we show 
how the asymptotic expansion for small correlation times can be obtained and we give some explicit 
results for simple noise processes. We also include a comparison between the asymptotic expansion and 
the exact numerical results to show the range of validity of the asymptotic expansion. Finally we give 
some conclusions. 

2 Definition of the model 

We consider the one-dimensional motion x(t) of a particle in a dissipative environment. The particle 
moves in a one-dimensional periodic potential V(x) with period L, and it is subject to a stochastic force 
z(t) that is correlated in time. In the over-damped regime, the motion of the particle is described by a 
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Langevin equation of the form. 

-=f( x ) + V2fW)+z(t). (2.1) 

We use units where the friction constant is unity. The first term on the right hand side is the force 
f(x) = — Qj- due to the potential V(x). The second term describes thermal fluctuations, t,(t) is a white 
noise with zero mean and (^(/)^(/')) = S(t — t'). In the later part of the paper we will discuss the zero 
temperature case, T = 0. The additional noisy force z(t) has zero mean, (z{t)) = 0. We assume that it is 
described by a Markov process with an infinitesimal generator M z . The probability density p(z,t) of this 
Markov process satisfies the Fokker-Planck equation 

^M=M^fc«). (2-2) 

We discuss a class of Markov processes for which the generator M z is described by its eigenvalues and by 
certain properties of its right eigenfunctions, namely 

M z § n (z) = -X s <Sfn(z), z§q{z) =Yo,i<I>i(z), 

zi(z)=yfl !fl +ii+i(z)+T«,»-ii-i(z) n = \,2,... (2.3) 
The eigenvalues — X„ obey 

h n < if n < m, Xq = (2.4) 

(j)o (z) is the stationary distribution of z. Due to the recursion relations the eigenfunctions (j)„ (z) can be 
written as (t>„(z) = g n (z)§o(z) where g n (z) are orthogonal polynomials with respect to the weight function 
<t>o(z), 

dzgn(z)g m (z)tyo(z) = 8„, m . (2.5) 

This class of Markov processes is very general. It contains many processes that occur in typical situa- 
tions such as the Ornstein-Uhlenbeck process, the dichotomous process, sums of dichotomous processes, 
and kangaroo processes. The correlation time x of the process is related to the smallest non-vanishing 
eigenvalue of the noise process via x = 1 . 

The joint probability density p(x,z,t) for the two stochastic variables x(t) and z,(t) obeys a Fokker- 
Planck equation of the form 

dp{X dt Z,t) =-y x (f(x)+z)p(x,z,t)+M z p(x,z,t) (2.6) 

We discuss only stationary properties of this equation, therefore the left hand side is put to zero and ( |2.6| ) 
becomes an equation for the stationary probability density p(x,z)- Due to the periodicity of V(x) we 
assume p(x,z) to be periodic in x with period L. To solve the stationary Fokker-Planck equation it is 
useful to expand p(x,z) in terms of the right eigenfunctions of M z . 

p(x,z) =poW^o(z) + £(-l)>„(z)^W. (2.7) 
This yields the recursion relations 

J = f(x)p (x)-y Q Ap\(x) (2.8) 
Yo.iPoM = hpi{x) +f(x)p[(x)-y 1 ,2p' 2 (x) (2.9) 



3 



Yn-\,nP' n -\{x) = KPn(x) + f(x)p'„ (x) - ln.n+\P n+ \ (x), U>\. 



(2.10) 



for the functions p n {x). As additional conditions we have the normalization of po(x) and the periodicity 

Of p n {x), 
L 

P0 (x)dx = 1 , p n {x)= p n (X + L). (2. 1 1) 



po(x) is the stationary distribution of x(t). It is one of the quantities we are interested in. A second and 



more important quantity is the integration constant J in (2.8), it is the stationary current. If the potential 
V(x) has inversion symmetry, the current vanishes. In the general situation, where V(x) has no inversion 
symmetry, one obtains generically a non- vanishing current. 



To solve the recursion relations ( |2.8| - |2.10| ) for a general Markov process M z approximations are nec- 
essary. One has to truncate the recursion at some large value N, i.e. one has to put Jn.n+i = 0. For some 
processes, e.g. for a sum of dichotomous processes, one has Jn,n+i = f° r some finite value of N and the 
approximation is not necessary. In principle it is then possible to solve the recursion relations numerically. 
Here, we restrict ourselves to a class of piecewise linear potentials. This is a standard assumption [g], 
such a potential is often called a ratchet-like potential. If the potential is piecewise linear, the force f(x) 
is piecewise constant. The differential equations ( pT8 - 2.1Q| ) are equations with constant coefficients that 



can be solved analytically. The remaining algebraic problem is to solve some continuity conditions for 
the functions p n (x) at the points where the force jumps from one constant value to another. In that way it 
is possible to express the stationary current as a ratio of two determinants. It can thus be calculated easily 
numerically for various parameters of the system or even analytically for small values of N. The general 
procedure has been described for finite temperatures in detail in One can obtain very accurate results 
in the whole parameter regime and for various noise processes. Thereby one observes that several ana- 
lytical approximation schemes do not work in different regimes of the parameter space. For instance, the 
standard expansion for small correlation times x breaks down. The physical reason is that for sufficiently 
small values of z the particle cannot escape one of the minima of the potential. At zero temperature, this is 
only possible if \z\ is larger than the largest value of |/(x)|. Since the process z(t) has a finite correlation 
time, z remains for a while in the region where \z\ < \f{x)\ and the particle cannot move apart. For a gen- 
eral potential V(x) this yields divergencies in the stationary distribution po(x), in the case of a piecewise 
linear potential, the stationary distribution contains contributions of the type Wi&(x — xi), where x\ are the 
positions of the minima of V(x) and w, are some weights. Such an effect occurs always if there is a finite 
probability to have a stochastic force z(t), that fulfils the condition /(jc; — 0+) < z < /(x, + 0+). It occurs 
especially in the case where Jn.n+i = for some finite N. Mathematically it turns out that the stationary 
distribution po(x) and the current J are nonanalytic functions of x in some regions of the parameter space. 
This is the reason why the x-expansion fails for some noise processes. 

In the next section we describe what modifications to the method in [Q] are necessary so that T = 
can be treated as well. The main goal is to obtain correct zero temperature results for the stationary 
behaviour of the model and to derive a correct asymptotic expansion for small x. 

3 The method at zero temperature 

In the rest of the paper we restrict ourselves to the simplest nontrivial case for the force f(x); we assume 
that f(x) takes to different values. To be precise, we let 

f(x)=f l if 0<*<Z4, 

f(x)=h if U<x<L. (3.1) 
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We let L2 = L — L\ . Due to the periodicity of V (x) we have f\L\+ ^2^2 = 0. We assume that V (x) has a 
minimum at x = 0, which means f\ < 0, fi > 0. Due to the discussion above, po(x) has the form 

po(x)=po(x)+W Q 5(x). (3.2) 

po(x) contains no further 6-contributions. In the same way we obtain for the functions p' n {x) 

p' n {x)=p' n (x)+W^{x). (3.3) 

The coefficients Wq and W' n can be related to each other using the differential equations 2.8- 2.1C| . One 
obtains 



2k-\ 

wl 



2k 





(-1) 



k Y2fc-2,2*-l 



Wo 

Yl,2 



J2k-\,2k 

For po(x) and p' n (x) we make the ansatz 



Po,i(x) 

p\A x ) 



r 



+ 1 



(3.4) 
(3.5) 



(3.6) 
(3.7) 
(3.8) 



The index i takes the two values i = 1,2, according to the two regions where f(x) = f\. Inserting this 
ansatz in the differential equations, we obtain a generalized eigenvalue problem to determine the coeffi- 

01 (r) (r) 

cients a- ,a ni , and b n • . This is the usual procedure for a system of coupled linear differential equations 
with constant coefficients and has been described for the present case in detail in We obtain using a 
vector notation 



/ b ,i \ 



b 2 j 
V : J 



1 ? ^ 

Jyo.i 

hf, 


V : J 



(3.9) 



(r) (r) 

The generalized eigenvalue problem for a- and a n ■ is 
where 



(3.10) 





( 





. 


.. \ 







*1 


. 





































B,- 



/ -fi Yo,i 
70,1 -fi Yi,2 

Yi,2 -fi 



(3.11) 



.(°) 



a) is zero, the other eigenvalues can be determined either be explicitely solving the eigenvalue problem 
or by using certain continued fractions as described in [Q]. To finally solve the problem, one has to 
determine the coefficients c r j, the weight Wo of the 6-contribution in po(x), and the current /. These are 
2N + 2 unknown variables. The continuity of p n (x) at x = L\ and certain jump conditions of p n {x) at 
x = involving the coefficients W' r for n > 1 yield 2N equations, in addition we have the normalization 
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of po(x). The last condition can be obtained from the continuity of the current density at a fixed value of 
z. For finite N we can write p(x,z) in the form 



p(x,z) = £ pM(x)&(z-Zk) ■ (3.12) 

k=Q 



for a process z(t) that takes N +1 different values Zk- P^ k '(x) can be related to po(x) and p' n {x), for 
instance 

jrpM(x)=p Q (x). (3.13) 
For each PW (x) we have a corresponding current 

./*(*) :=**(*)/*(*)■ (3-14) 

We introduced /*(jc) := /; + z\ if jc lies in the region where f(x) = f t . The current J k (x) has to be 
continuous at x = L\. If \zk\ < min,(|/;|) and sign(zt) = — sign(/|) this yields J k (L\) = 0. With this 
equation we have found the last condition we need to determine the 2N + 2 unknown variables. In this 
way the calculation of the current J and the stationary distribution po(x) has been reduced to the purely 
algebraic problem of solving 2N + 2 linear equations for 2N + 2 variables. As in [Q], the current 7 can be 
expressed in closed form using determinants. 



4 Results 

The results presented in this section have been obtained by solving the algebraic problem mentioned 
above numerically. We present results for two classes of processes, namely for a sum of dichotomous 
processes and for some kangaroo processes. 

4.1 Sums of dichotomous processes 

A dichotomous process z{t) is a process where z(t) takes two values, ±zo- Summing up Af such processes, 
we obtain a Markov process with a stationary distribution of the form 



Hz) = ^ E , )Hz-(N-2i)zo) ■ (4.1) 




The parameters that characterize the process are 

Xn = -n/%, n = 0,...,N, (4.2) 



y n ,n+i = V(n + l)(N-n)zo, n = 0,...,N. (4.3) 



In order to have a value for (z 2 ) = Jq l that is independent of N, we choose zo = y/ VN where y = yjD/x. 
D is the noise strength of the process. In the limit N — > °°, the sum of dichotomous processes yields the 
Ornstein-Uhlenbeck process [0]. 



The eigenvalue problem ( |3. 10| ) can be solved explicitely for a sum of dichotomous processes. This 



has been shown in The eigenvalues a\ k ^ are 

a ; (1) = 4^4 (4.4) 
' _ ff 
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a? kU2k+l) = ^ ? ,for2k+KN (4.5) 

2 (N -2k) 2 f-Nff ' 

and in addition 

aj --ttf (4.6) 

for even A". The coefficients 7, Wo, and c r j can be determined numerically as described above. Let us 

first discuss the results for small correlation times x. Since y 2 = D/%, X = 1/x, the eigenvalues CC-^ have 
singularities oc x~ 1//2 and oc x _1 . The second type of singularity is present only if N is even. Furthermore, 
for small x and N odd, \z(t)\ is always larger than \f(x)\ and we expect that Wo = 0, whereas for even 
N, z(t) may be zero and therefore we expect Wo > 0. The behaviour should be very different depending 
on whether N is even or odd. First we will show some results for the current. In Fig. la the current is 
plotted as a function of x for N odd. The current is always positive. The corresponding result for even 
N is shown in Fig. lb. The current is negative for small x and changes the sign when x becomes larger. 
This has already been observed for very low temperatures in [Q]. We will come back to this point in the 
next section, where we discuss an expansion for small x. Let us mention that a negative current for even 
N and small x occurs only if the noise strength D is sufficiently large. For larger values of x, the current 
is shown for some processes in Fig. 2. The current has cusp-like maxima that occur near such values of 
X, for which one of the discrete values of \z(t) \ becomes smaller than one of the discrete values of |/(x)|. 
The exact position of these values of x is indicated in the figure. It is interesting to see that the maxima 
don't lie exactly at this value of x but slightly below. With increasing number of maxima, i.e. for higher 
values of x this effect increases as well. Up to now we have no physical interpretation for it. 

Whenever x becomes larger than one of these values of x, the weight Wo of the 8-peak in the stationary 
distribution po(x) jumps discontinuously to a higher value. This is shown in Fig. 3 for the same parameter 
as used in Fig. 2. The weight Wo shows the expected behaviour and is unity for large x. The coefficients 
c r j can be calculated as well, they enter in the expression for the stationary distribution po(x). Typical 
results are shown for various processes in Fig. 4. x is very small, so that a 8-contribution (which is not 
shown) occurs only for even N. One observes that for both cases, even or odd values of N, the stationary 
distribution po(x) has discontinuities. We will come back to that point later. The minimum at x = 0.8 
is located at the maximum of the potential. It is strong for N = 2 and becomes weaker when N is even 
and increases. The opposite behaviour is found for odd N, here the minimum is only weak for N = 1 
and becomes more and more pronounced for larger N. Similarly, the maximum at x = increases with 
increasing N if N is odd. 



4.2 Kangaroo processes 

Kangaroo processes are processes with X n = — 1 /x for all n > 0. They can be completely characterized 
by the stationary distribution p s t{z)- In this section we discuss results for kangaroo processes with a 
stationary distribution given in Q4.1[ ). According to the notations there, we call these processes K(N). This 
allows a direct comparison with the results for the sums of dichotomous processes. The only difference 
between the kangaroo processes and the sums of dichotomous processes is that for the latter case, only 
jumps of a magnitude ±zo in z(t) occur, whereas a kangaroo process has no restriction of the jumps of 
z(t). Since the stationary distributions are the same, the parameters y, un +i are the same. But unfortunately, 

(r) 

one cannot derive an analytic expression for the eigenvalues a- for N > 4. They have to be determined 
numerically, which is easily done. 

We show some of the numerical results for the current as a function of the correlation time x in Fig. 5a 
for N odd and Fig. 5b for N even. The inset in Fig. 5a shows that the current is negative for very small 
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X. This is in contrast to the sums of dichotomous processes. Furthermore, the region of negative current 
becomes larger, when D becomes smaller. For even N the current reversal is observed as for the sum of 
dichotomous processes, but the region of negative current is larger. Furthermore the current is always 
negative for small x and N even, whereas for the sum of dichotomous processes a negative current occurs 
only if the noise strength is large. Thus, the behaviour of the current as a function of x and D differs 
significantly from that one of the sums of dichotomous processes, although the stationary distribution is 
the same in both cases. Let us mention that the mechanism that produces a negative current for small 
X and N odd is not clear. Doering et al. ^ proposed that a negative current occurs due to a parameter 
called flatness, which is a property of the stationary distribution. Their argument cannot be applied in 
the present case, since the stationary distributions and therefore the flatness is the same for the kangaroo 
processes and the sums of dichotomous processes, whereas the sign of the current for small x is different. 



5 An expansion for small x 

The usual formal expansion for small x has already been discussed by Doering et al It can be ob- 
tained using a standard perturbational treatment [0] to solve the stationary Fokker-Planck equation. The 
unperturbed part is M z , which is of order x . The perturbation contains a term proportional to z„ it is 
thus x -1 / 2 . But since it contains only off-diagonal matrix elements in the basis in which M z is diag- 
onal (there are no parameters y„,„), one obtains finally an expansion in powers of x for the stationary 
distribution and for the current J. Since the current vanishes in the white noise limit, the first term in the 
expansion for J vanishes and / « x. An alternative way to obtain the same expansion is to use a formal 



operator continued fraction that can be obtained from the recursion relations (2.8-2.10). This has been 
shown in the appendix of [Q, where a general expression for the first order has been given and higher 
orders can be obtained straight forward. In our special case of a saw-tooth potential one obtains 

7 = x7i+x 2 7 2 + ... (5.1) 
with 




i 



eD - 1 ) I e d 




(5.3) 



This result generalizes the result given by Doering et al. [|3j]. It can be compared with the exact results 
of the sums of dichotomous processes and for the kangaroo processes shown in the last section. One 
observes that the formal expansion is correct for a sum of an odd number of dichotomous processes, 
if the correlation time x is sufficiently small. For a sum of an even number of dichotomous processes 
the result is not even qualitatively correct, the current has the wrong sign. Similar results hold for the 
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kangaroo processes. If N is odd, the expansion is correct for small x, whereas it is wrong for even N. The 
reason is that the stationary distribution po(x) contains a contribution Wq8(x) for even N and arbitrary 
small x, but not for odd N. This leads to non-analytic contributions to Pq(x) and J as functions of X if N 
is even, which are not taken into account in the perturbative x-expansion. This expansion is well defined 
only if po(x) is a differentiable function of x, which is true for odd N, but not for even N. 

A second possibility to obtain an asymptotic x-expansion for J and p$ (x) is to start from the set of 
linear equations that determine the variables J, Wo, and c r ,. The coefficients in this set of linear equations 
contain singularities of the type 

c 

e V t (type 1) 

c 

e (type 2) 

c 

e X (type 3) 

(k) 

where c is some constant. These singularities occur due to the behaviour of the eigenvalues a) for small 



X, see e.g. ( |4.4|416| ). For small x it is possible to divide those equation, which contain singularities of type 
1 by exp(+c/^/x). We can then neglect all the terms that vanish faster than any power of y/x. For small N 
it is possible to solve the remaining set of linear equations analytically. For a single dichotomous process 
this yields 

jO)(x) =7 1 (1) x + 7^ 1) x 2 + 0(x 3 ) (5.4) 
where 

/« = fi-fi (55) 

&(j- i -j-\ 2 (eT>-\\(elT 



Ti 

AD _ g + fi 



f 4 _l f 4 

l l V / i \ 2 / i \ 2 



I~h leT '~ u [l '~ T '- ] 



i (ft+f?) 



fx -ft 



\c--\ \ [C 77-1 



1 l^ 2 



(5.6) 



i 



As expected, this result agrees with the formal x-expansion described above, which was correct for a 
single dichotomous process. For a sum of two dichotomous processes we obtain 

J {2) = l -4 l \l-D)x + 0(z 3 ). (5.7) 

This expression differs from the formal x-expansion, but it agrees well with the exact numerical results 
for small x. The current is negative for small x and D > 1. For N = 3 one would again expect a result that 
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agrees with the formal x-expansion, but this is not true. We obtain 




(5.8) 



(5.9) 



2^6 ( e h _ i) (el _ i) 



The additional term <=c \i was not present in the formal x-expansion. It is responsible for the very small 
region of validity of the formal x-expansion for N = 3. Similar terms occur for any odd N. The reason 
is that the stationary distribution po(x) is not a differentiable function of x. It has discontinuities at the 
extrema of the potential, see Fig. 4. Therefore the perturbative x-expansion is not well defined. 

Finally one can look at the kangaroo process K(2), for which the formal x-expansion was incorrect 
too. We obtain 



which is in agreement with the exact numerical results. A comparison between the numerical results and 
the asymptotic expansions is shown in Fig. 6. 

6 Conclusions 

We discussed in detail the behaviour of the static properties of one-dimensional models for noise induced 
transport at zero temperature. For our discussion we used a sawtooth potential, but the method and 
the results can be generalized easily to any piecewise linear potential. If the support of the stationary 
distribution of the stochastic force z(t) is finite, one observes non-analytic contributions in the induced 
current and in the stationary distribution of the coordinate of the particle. In a piecewise linear potential, 
the stationary distribution of the coordinate po(x) contains contributions of the form w;8(x — x,). Such 
contributions occur for sufficiently large values of x, for some processes even for all x > 0. If the (j)o (z) 
and the force f(x) = —dV/dx are continuous, one can not expect a 8-contribution to p(x). But even then 
one expects singularities in p(x). Let / m ; n and / max be the minimum and the maximum of f(x), and let 
us assume without loss of generality that |/ mm | < |/ max |- Let us now take a stochastic force z(t) that takes 
only a finite, discrete set of values zi- If < Zi < —fmin, the particle moves to the right until z, = f(x). 
At these special values of x singularities in p(x) may occur, depending on the behaviour of f(x) near 
such a point. A similar situation occurs if > Zi > —/max- In the more general case where the stochastic 
force can take all values within an interval [— zo,Zo], similar singularities may occur, again depending on 
the form of f{x). Therefore the results we presented are relevant for a large class of situations. Let us 
mention that one should not expect such singularities in the case of a Gaussian noise process, since in that 
case the stochastic force may be arbitrarily large. 

A consequence of the non-analytic behaviour of p(x) is that the usual perturbative x-expansion breaks 
down. In the n-th order of this expansion a derivative of order (n + 2) of p(x) occurs. If this derivative 
doesn't exist, it is clear that the x-expansion is not defined. In contrary, the asymptotic x-expansion we 
derived for some special cases is always well defined. 



jK(2) = J m x+0{z 3 ) 

Jf® = -Dj[ l \ 



(5.10) 



(5.11) 
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Figure captions 

Fig. 1. The current / as a function of the correlation time X for a sum of ./V dichotomous processes, a) N 
odd and D = 10 AV and Li/L 2 = 4. 
b) The same as in a) for N even. 

Fig. 2. The stationary current for the sum of N = 1,2,3 dichotomous processes as a function x on a larger 
range of x. The tags indicate the correlation times, for which is kzo = fi- The curves for 
N = 1,2,3 have maxima near X\k,^2k^3k- The parameters are the same as in Fig. 1. 

Fig. 3. The weight of the 8 distribution as a function of X for the sum of N dichotomous processes for the 
same parameters as in Fig. 1. 

Fig. 4. The stationary probability distribution of the position coordinate in a normed period and with 
L\/Li = 4 for the sum of N dichotomous processes. N = \ (solid line), N = 2 (dashed line) N = 3 
(short dashed line) and N = 4 (long dashed line). The 5-contribution is not shown. 

Fig. 5. a) The current J as a function of the correlation time x for the kangaroo processes K(N), N, odd 
with the same parameters as in Fig. 1. The inset shows the same plot on a finer scale, 
b) The same as in a) for N even. 

Fig. 6. Comparison of the linear approximation (solid lines) of the current as a function of the correlation 
time with the exact results (dashed lines) for the sum of N dichotomous processes and the kangaroo 
process K (2). 



Fig. l.a) 



y/i<r 2 Av 



10 









N=3 






N=5 














N=9 


N=1 / 










N 


=7 





7/AV 



Fig. l.b) 



0.5 



0.3 



0.1 



-0.1 



-0.3 



N=8 










_ — N=4 




N=2 







0.0 



0.1 

t/AV- 1 



0.2 



Fig. 2 




14 



Fig. 3 




15 



Fig. 4 




16 



Fig. 5a) 



7/10 3 AV 




Fig. 5b) 



J/AV 



N=10 




0.0 



0.1 



0.1 



0.2 



t/AV- 1 



17 



Fig. 6 




18 



